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ABSTRACT 

We report about very high energy (VHE) 7-ray observations of the Crab Nebula with the MAGIC telescope. 
The 7-ray flux from the nebula was measured between 60 GeV and 9 Te V. The energy spectrum can be de- 
scribed by a curved power law dF /dE = fo (is/SOOGeV) [ fl+Mo Sio( £ / 300GeV )] w j m a g ux norma ]i za tion /o of 
(6.0 ± 0.2 stat ) x lO-^cirrV'TeV- 1 , a = -2.31 ± 0.06 stat and b = -0.26 ± 0.07 stat . The peak in the spectral en- 
ergy distribution is estimated at 77 ± 35 GeV. Within the observation time and the experimental resolution of 
the telescope, the 7-ray emission is steady and pointlike. The emission's center of gravity coincides with the 
position of the pulsar. Pulsed 7-ray emission from the pulsar could not be detected. We constrain the cutoff 
energy of the pulsed spectrum to be less than 27 GeV, assuming that the differential energy spectrum has an 
exponential cutoff. For a super-exponential shape, the cutoff energy can be as high as 60 GeV. 
Subject headings: acceleration of particles — gamma rays: observations — pulsars: individual (PSR 
B0531+21) — radiation mechanisms: non-thermal 
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1. INTRODUCTION 

The Crab Nebula is the rem nant of a supernova explosion 
that occurred in AD 1054 (e.g. lCollins et al.|[l999l and refer- 
ences therein) at a distance of ~ 2kpc. It is one of the best- 
studied non-thermal celestial objects in almost all wavelength 
bands of the electromagnetic spectrum from 10~ 5 eV (radio) to 
nearly 10 14 eV (7-rays). There is little doubt that the engine 
of the nebula is the pulsar PSR B053 1+2 1 (hereafter Crab pul- 
sar) in its center. 

In very high energy (VHE) 7-ray astronomy the Crab Neb- 
ula was first detected with large significanc e at TeV en- 
ergies by the pioneering Whipple telescope (Week eset al.l 
1989). The nebula turned out to be the strongest source 
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of steady VHE 7-ray emission in the Galaxy. It is there- 
fore used as the standard "calibration candle" for ground- 
based 7-ray experiments. Apart from testing the perfor- 
mance of 7-ray instruments another aim of measuring the 
Crab Nebula is to increase the measurement precision of 
the Crab Nebula flux and the energy range covered. These 
continuous efforts provide insights necessary for the under- 
standing of the very details of the emission mechanisms 
of VHE 7-rays in the Crab Nebula and pulsar. Impor- 
tant questions remain to be answered concerning the emis- 
sion mechanisms of the nebula, of the pulsar at GeV ener- 
gies, and of the nebula around PeV energies. Since its dis- 
covery, detailed studies of the VHE emission energy spec- 
trum, ranging from s everal hundred GeV up to 80 TeV, have 
been carried out (e . g.lAkerlof et al.lll990HVac anti et al J| 199 l b 
Baillon et alj|1991t iGoret et al.lll993t iKonopelko et alJll996t 
Tanimori et al. 1998c iHillas et alj 119981: lA menomo ri et al.l 
1999t iMaiumdar et al.1 12002c lAharonian et alj 12004 120061) . 
Between 10 and ~ 200GeV, observations are sparse. A few 
results are provided by converted solar co ncentrator arrays 
that use the wave front s ampling technique (jOser et alj|2001b 
de Naurois etail l2002h lArqueros etail [2002; Smit h et all 
2006). However, the wave front sampling suffers from rel- 
atively large uncertainties in the calculation of the effective 
area and of the energy, as well as from poor 7/hadron discrim- 
ination, which make it difficult to perform differential flux 
measurements. 

A good explanation of the nebular dynamics and the ob- 
served energy spectrum below GeV energies can be ob- 
tained with the magne t o-hyd rodynamic model suggested 
first by [Re es and Gunnl d 19741) and developed further by 
iKennel and Coronitil d!984albl) . In this framework, the pul- 
sar provides a continuous flow of charged particles (pulsar 
wind) with Lorentz factors of 10 — 10 7 . A standing reverse 
shock forms where the wind ram pressure balances the total 
pressure of the nebula. Wind particles accelerate in the shock 
to ultra-relativistic energies and subsequently lose their en- 
ergy by synchrotron emission. The presence of synchrotron 
emission up to a few hundred MeV in conjunction with the 
observed 7-ray spectrum at TeV energies shows that particle 
acceleration takes place up to energies of ~ 10 15 - 10 16 eV. 
From the total luminosity of the synchrotron emission, it can 
be inferred that about 10% of the pulsar's energy loss rate is 
converted into the kinetic energy of particles. 

Above 1 GeV the dominant source of 7-ray emission is 
most likely inverse Compton (IC) scattering of synchrotron 
photons by the synchrotron-emitting electrons in the s hocked 
wind region (synchrotron s erf-Compton model fSSC|: | Gouldl 
fl965t IWeekes et alJ fl989t Ide Jager and Hardind [T992Ik To 
explain the observed VHE 7-ray spectrum, several other 
seed photon fields are also believed to contribute to the in- 
verse Compton scattering, namely, far-infrared excess, cos- 
mic microwave bac kground, and millimeter-photons (e.g. 
lAharonian et~aill2004l) . 

Although the IC-mechanism gives a good description of 
the observed energy spectrum between 500 GeV and about 
10 TeV, other processes may contribute in part to the VHE 
7-ray emission. It is likely that a significant fraction of 
the mechanical energy lost by the pulsar is taken away by 
a hadronic component in the wind. Following interactions 
of this component with the interstellar medium, VHE 7- 
rays are emitted by decaying 7r°s, whi ch modify the energy 
spectrum at TeV energies and beyon d ( Atovan and Aharonian 
1 19961 iBednarek and Profheroelll997[ fBednarek and Bartosik 



120031: lAmato et alJ[2003h . lAtovan and Aharonianl d 19961) dis- 
cuss the possibility of an "amplified" bremsstrahlung flux at 
GeV-energies, which could account for the discrepancy be- 
tween the measured GeV 7-ray flux and predictions within 
the SSC-framework (de Jage r et all 119961) . If this is true, 
one should observe, in good approximation, a power-law 
spectrum between 100 GeV and 10 TeV with a spectral in- 
dex 2.5-2.7. Another mechanism to be mentioned is IC- 
scattering of relativistic electrons in the unshocked pulsar 
wind. If the target photons are emitted by the pulsar, a 
pulsed compo nent could extend to 7-ray en ergies of sev- 
eral 100 GeV (Bogovalov and Aharonian 2000). An indepen- 
dent measurement in the intervening region between 60 and 
400 GeV would constrain further the parameters of various 
models. The MAGIC imaging atmospheric Cerenkov tele- 
scope (IACT) has a low-energy trigger threshold (^50 GeV) 
and is currently the only experiment capable of exploring this 
energy regime. 

The spatially resolved morphology of the nebula is of a 
complex nature. Its size in optical bandwidths is about 4 x 
6 . Due to synchrotron losses, the high-energy electrons will 
have shorter cooling times and only the lower energy elec- 
trons will reach out farther into the nebula. Thus, the effective 
source size is expected to shrink with increasing energy of 
the radiation. The radio emission is expected to extend up 
to and beyond the filaments optically visible, whereas X-ray 
and multi-TeV 7-rays should be produced in the vicinity of 
the shock. On the other hand, the expected source size would 
increase if the presence of an ionic component is established. 
In a special study the HEGRA collaboration concluded that 
the rms size of the V HE 7-ray emission reg ion is < 1.5' for 
energies above 1 TeV (Aharon ian et al.l20 00). A similar study 
at energies below 1 TeV has not been performed up to now. 

The Crab pulsar is a source of pulsed radiation and 
has been detected up to GeV energies. The Crab pulsar 
has a period of 33 ms and was first discovered in radio 
dStaelin and Reifensteinl [l968) and shortly afterward s as the 
first pulsar in the optical domain dCocke et alj[l9 69). Since 
then, pulsed emission from the Crab pulsar has been de- 
tected at all accessible energies up t o 7-rays (for a compila - 
tion of the broad band emission see Thom pson et al.l 1999). 
EGRET detected the Crab pulsar in 7-rays up to energies 
of 10 GeV with a hint of a cutoff in the energy spe ctrum 
in the highest energy bin at 6GeV dNolan et alJ ll993). De- 
spite various efforts, observations from the ground at higher 
energies have, so far, failed to detect pulsed emission (e.g . 
MusauerdfT999t lAharonian et al.lll999t iLessard et afl 120001: 
de Naurois et al . 2002; A haronian et alJl2004|). Some experi- 
ments reported episodic pulsed emissi on (Gibson etai1 [l982: 
iBhat et al.lll986h lAcharva etail 1 19921) and persistent p ulsed 
emission over a 1 vr period dDownthwaite etail [19841) . but 
these observations have not been confirmed by other experi- 
ments. 

The high-energy emission from the pulsar is assumed to be 
due to curvature and synchrotron radiation from relativistic 
charged particles that are forced to move along magnetic field 
lines inside the magnetosphere of the pulsar. The question 
of where the particles are being accelerated is the subject of 
ongoing theoretical activities. In the two most popular mod- 
els, the production of electrons and positrons and their accel- 
eration tak e place either above the polar cap of the neutron 
star (e.g. lHarding et al.lll978t iDaughertv and Hardindll982l) 
or in outer gaps in between the null surface and the light 
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cylinder of the magnetos phere (e.g. ICheng et ail [l986b a; 
IChiang and Romanilll992l) . We should not omit the slot-gap 
model, which pl aces the acceleration zone at the outer rim 
of the polar cap dArons| [T983; Muslimo v and Harding||2003l) . 
These models differ in the predicted shape and cutoff of the 
energy spectrum at the highest energies. A measurement of 
the turnover in the spectrum would shed light on the possible 
sites of the particle acceleration and can constrain models. 

In this paper we report about the observation of the Crab 
Nebula and pulsar with the MAGIC telescope between 2005 
October and December. After describing the MAGIC tele- 
scope and the performed observations in §§|2] and [3] we 
present the analysis chain in $4] In §|5] we present our re- 
sults on the steady emission comprising the differential flux 
between 60GeV and 9TeV, a study of the morphology of 
the emission regio n, and the measured integral flux above 
150GeV. In jj5.4l we present the results of our search for 
pulsed emission, and we close with a discussion of our results 
ing6] 

2. THE MAGIC TELESCOPE 

The MAGIC (Major Atmospher ic Gamma Imaging 
Cerenkov) telescope (see Lorenz 2004) is located on the Ca- 
nary Island, La Palma (2200 m asl, 28.45° north, 17.54° 
west). MAGIC is currently the largest single-dish IACT 
It has a 17 m diameter tessellated reflector consisting of 
964 (0.5 x 0.5) m 2 diamond-milled aluminium mirrors, which 
are grouped onto support panels in units of four. Depend- 
ing on the elevation of the telescope, the position of every 
panel is adjusted by computer-controlled actuators of the so- 
called Active Mirror Control, thus providing optimal focus- 
ing. About 80% of the light from a point source is focussed 
within a radius of 0.05° in the focal plane. The MAGIC tele- 
scope is focused to a distance of 10 km — the most likely loca- 
tion of the shower maximum for 50 GeV 7-ray air showers at 
small zenith angles. 

The faint Cerenkov light flashes produced in air show- 
ers are recorded by a camera comprising 577 photomulti- 
plier tubes (PMTs). The inner part of the camera (radius 
- 1.1°) is equipped with 397 PMTs (type ET 9116A; 1" 
from Electron Tubes Ltd [ET]) with a diameter of 0.1° each 
. The outer part of the camera is equipped with 180 PMTs 
(type ET 9117A; 1.5" 0) of a diameter of 2°. The cen- 
tral P MT is modified for optical pulsar studies dLucarelli et alJ 
120051) . Hollow hexagonal non-imaging light concentrators 
(often called light catchers) are placed in front of all photo- 
multipliers to compensate for the dead space between them 
and in order to shrink the observation solid angle to the reflec- 
tor. The entrance window of the PMTs is coate d with a diffuse 
lacqu er doped with a wavelength shifter (WLS; Paneq ue~et al.l 
2003). The combination of the hemispherically shaped PMT, 
the light catcher, the diffuse coating and the WLS results in 
a 15%-25% higher quantum efficiency compared to flat win- 
dow PMTs. For protection purposes (humidity, dust) a thin 
entrance window made of Plexiglas (type UG-218, with a UV 
cutoff around 290 nm) is placed in front of the camera. 

The PMTs have 6 dynodes and operate at a gain of roughly 
30,000 to slow down aging and damage from high currents by 
light during observations close to the Galactic plane and dur- 
ing moonshine. After amplification, the fast analogue signals 
are converted to optical signals and transported by 160 m long 
optical fibers to the counting house. There the optical signal 
is converted back and split. Part of the signal is routed to the 
trigger. The current configuration of the MAGIC camera has 



a trigger region of 2.0° in diameter dCortina et al.ll2005t) . This 
provides a 7-ray trigger collection area of the order of 10 5 m 2 
(at 200 GeV for a source close to zenith). Presently, the trig- 
ger energy range spans from 50-60 GeV (peak in the differen- 
tial trigger rates at small zenith angles for a 7-ray source with 
a spectral slope of -2.6) up to tens of TeV An event is trig- 
gered if the signals in each of 4 neighbouring pixels exceed 
a threshold of ~ 7photoelectrons (phe) within a coincidence 
time window of 6 ns. 

Before being digitized by an 8-bit, 300 MSamples/s FADC 
system, each signal is stretched to an FWHM of about 6 ns. 
The FADC continuously writes the digitized amplitude infor- 
mation into a ring buffer. In case of a trigger the digitization 
stops and the corresponding part of the ring buffer is written 
onto a disk. The dead time introduced by the readout is 25 /is. 
In order to expand the dynamic range to ~ 1000, the signal 
of every PMT is split into two branches, differing by a factor 
of 10 in gain. The higher gain branch is read out for a 50 ns 
time interval. When the signal amplitude exceeds a preset 
threshold, the delayed lower gain is routed to the same FADC 
channel and recorded in the following 50ns. Otherwise, the 
signal from the high gain branch continues to be recorded and 
is used to determine the pedestal offset of each PMT channel. 

The accuracy in reconstructing the direction of incoming 7- 
rays on an event-by-event basis, hereafter 7-ray point spread 
function or 7-PSF, is about 0.1°, depending on the energy. 
With the information provided by a starguider camera, mis- 
pointing is corrected to an absolute precision of about 1'. A 
7-ray source with an absolute intensity of ~ 2% of the Crab 
Nebula and similar spectrum can be detected with MAGIC 
within 50 hr at energies > 200 GeV on a significance level of 

5(7. 

3. OBSERVATIONS AND DATA SELECTION 

Observations of the Crab Nebula with MAGIC are con- 
ducted on a regular basis, as a means to monitor the perfor- 
mance of the telescope. In this report we restrict ourselves to 
the analysis of data obtained in the first observation cycle of 
the MAGIC telescope between 2005 October and December. 
The observations were performed in the so-called ON/OFF 
mode. The telescope was pointed towards the Crab pulsar 
(ON) for about 16 hr. An OFF source position, a sky region 
where no 7-ray source is known, was observed in the same 
range of zenith angles as the ON source. For the background 
estimation we used OFF data collected for over 19 hr. 

One of the main objectives of this analysis was to explore 
the lower energy range of accessible 7-ray energies. The en- 
ergy threshold of IACTs depends strongly on the zenith angle 
of observation. Restriction to events with low zenith angles 
provides the lowest possible energy threshold. Therefore, we 
select events with zenith angles < 20°. For any given night the 
data affected by technical problems or fluctuations in the data 
rate in excess of 10% were rejected. The atmospheric con- 
ditions were judged from the nightly averaged and publicly 
available atmospheric extinction coefficients from the nearby 
Carlsberg Meridian telescope 23 . Within the selected nights, 
the atmospheric light transmission changed by less than 5%. 
The nights with data that survived all the selection criteria, to- 
gether with the corresponding observation times, trigger rates, 
and zenith angle range, are listed in TableQ] The selected data 
sample comprises 14 nights amounting to a total ON observa- 
tion time of 955 minutes (~ 16 hr). 

23 See http://www.ast.cam.ac. uk~dwe/SRF/camc_extinction.html 
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TABLE 1 
Data Selected for Analysis 



Date Rates On Time Zd-Range 

(MJD) (Hz) (minutes) (deg) 

53,648 130 73 7-23 

53,655 115 100 7-19 

53.671 122 105 7-23 

53.672 105 61 8-20 

53,679 115 51 7-20 

53,684 95 50 7-20 

53,707 98 53 7-11 

53,709 105 48 7-10 

53,711 105 48 7-10 

53,713 108 50 7-14 

53,727 92 44 7-13 

53,729 97 61 7-13 

53,731 100 107 7-14 

53,735 92 104 7-15 



NOTE. — The second column lists the average event rate after tail cuts and 
a cut in SIZE> lOOphe (for the definition of SIZE see text). 

4. DATA ANALYSIS 

The data analysis was carried out using the stan- 
dard MAGIC analysis a nd reconstruction software MARS 
dBretz and Wagner 2003). After removing faulty and unstable 
camera channels, which amount to 3%-5% of the total num- 
ber of PMTs, the signal amplitudes, extracted with the digital 
filtering method dAlbert et alj|2006|). were converted to ph e 
by using the F-factor method (Mirzov an and Lorenzlll997l) . 
Using calibration events recorded interlaced to normal events 
(Gaug et al. 2005), the conversion factors were updated every 
10 s. There is a 10% systematic uncertainty in the calibration 
that directly propagates to the uncertainty of the event energy 
scale (point 7 in Table |9]below). Time offsets between pixels 
are corrected with a precision of better than 1 ns. 

After calibrating the data, pixels with faulty reconstructed 
signal amplitudes and times were rejected and the correspond- 
ing amplitude and time information was interpolated from 
the signals and times of neighbouring pixels. Before im- 
age parametrization, a tail-cut cleaning of the image was 
performed, requiring signals higher than a pre-defined ab- 
solute amplitude level and time coincidences (3.3 ns) with 
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FIG. 1 . — Parameterization of a shower image with Hillas-parameters. 
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FIG. 2. — Energy distribution of MC-7-ray events with SIZE> 100 phe for 
a simulated 7-ray source with spectral slope -2.6. 

neighbouring channels. The time coincidence effectively sup- 
presses pixels containing only a signal from the night sky. 
For most of the analysis the minimum required pixel con- 
tent is 6 phe for so-called core pixels and 4 phe for bound- 
ary pixels. For the morphology studies the minimum pixel 
contents were raised to 10 phe (core) and 8 phe (boundary), 
respectively, which improves the angular resolution, albeit at 
the expense of an increased analysis threshold. 

Every cleaned event was parameterized by a principal com- 
pone nt analysis, commonly referred to as Hillas parameteriza- 
tion dtjillas 1985). The parametrization was later used to sep- 
arate between 7-ray event candidates and background event 
candidates. The Hillas-parameters DIST, LENGTH, WIDTH 
and also ALPHA are illustrated for a recorded shower image 
in FigureQ] Another useful parameter is the SIZE of a shower 
image, the intensity of the image after image cleaning in units 
of recorded photoelectrons. Note that SIZE depends on the 
applied tail cuts. SIZE is a good estimate of the primary 
particle energy, provided that the shower impact distance to 
the telescope principal axis is below ~ 120 m. An event pre- 
selection was performed by discarding event candidates af- 
fected by noise and pick-up (e.g. car flashes) and event can- 
didates with a low number of pixels (after tail-cuts typically a 
minimum number of 5 core pixels were requested). In addi- 
tion, an image-SIZE of at least ~ 100 phe was requested. Fig- 
ure |2] shows the energy distribution of simulated 7-ray events 
with SIZE > lOOphe. The distribution peaks at an energy of 
75 GeV. The simulated 7-ray source has a powerlaw spectrum 
with an index of -2.6. 

A Monte Carlo (MC) simulation properly describing data is 
a necessary requisite for a ground-based 7-ray experiment. In 
Figure[3]the image parameter WIDTH is shown for simulated 
7-rays and 7-rays extracted from data. The four panels are for 
consecutive bins in SIZE covering the entire range of analyzed 
7-ray energies. For an unbiased comparison, loose cuts have 
been applied in the 7/hadron-separation (explained in the next 
section). The 7-ray excess was obtained by subtracting the 
scaled distribution of the OFF-data sample from the distribu- 
tion of the ON-data sample. The scaling factor was found by 
normalizing the |ALPHA|-distributions of both samples be- 
tween I ALPHA I =30° and 70°. The comparison was done by 
selecting events with small |ALPHA| (typically less than 10°). 
A small |ALPHA| is expected for 7-rays from the Crab Neb- 
ula, as explained later. The agreement of the MC simulated 
distributions and the distributions extracted from data is ac- 
ceptable in all four SIZE-bins. However, it is evident from the 
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figures that the agreement worsens at large SIZE. In the anal- 
ysis, a possible bias introduced by this behavior is avoided by 
applying rather loose cuts at energies above several TeV. 

4. 1 . Gamma Hadron Separation 

Only a small fraction between 10~ 3 and 10~ 4 of the recorded 
data are 7-ray showers. The major fraction of the recorded 
events are cosmic -rays of hadronic origin. This unwanted 
background has to be suppressed offline. For this purp ose we 
applied a multivariate method, the random forest (iBreimanl 
1200 It iBock etal.ir2004t lAlbert et al.ll2007bl) . which uses the 
image-parameters to compute the HADRONNESS of an 
event. The HADRONNESS of an event quantifies the prob- 
ability of an event to be "7-like" or "hadron-like." The ran- 
dom forest is trained with MC simulated 7-ray events and ei- 
ther with simulated hadronic cosmic-ray events or, as in the 
present case, with background events recorded by MAGIC. 
In this study, we used for the training of the random forest 
the image parameters SIZE, DIST, WIDTH and LENGTH, as 
well as the third moment along the major axis of a shower im- 
age, and a parameter describing the CONCENTRATION of a 
shower image. 

Figure g] shows the HADRONNESS for MC simulated 7- 
ray showers (red) and for recorded background (blue) as a 
function of SIZE. A clear separation between both popula- 
tions is visible for SIZE > 300 phe, corresponding to 7-ray 
energies > 150GeV. Below 300 phe both populations start to 
overlap until, at ~ 200phe (~ lOOGeV), no more separation is 
possible. The bottom panel in the figure shows the maximum 
quality factor Q (e-y/y/ei) for each SIZE interval obtained for 
an optimized HADRONNESS cut. Parameter e 7 is the frac- 
tion of retained 7-ray events and cb the fraction of retained 
background events. The corresponding HADRONNESS cut 
is shown by the stars in the top panel of the figure. 

The random for est method is also u sed to estimate the en- 
ergy of each event (lAlbert et al.ll2007bl) . An energy resolution 
of ~ 25% is achieved for events with energies > 200 GeV. The 
energy resolution reduces to ~ 40% around 70 GeV. 

4.2. Signal Extraction 

After 7/hadron-separation and energy estimation, the 7-ray 
signal is extracted from an | ALPHA | -distribution. ALPHA 
is the angle between the major axis of the recorded shower 
and the vector connecting its center of gravity (CoG) with the 
source position in the camera plane (cf. Figure [T|l. Shower 
images of 7-rays from the source point with their major axes 
toward the source position in the camera and appear as an ex- 
cess at small values in the | ALPHA | -distribution. Figure [5] 
shows two I ALPHA| -distributions from the data (black: ON- 
data; red: OFF-data). The left panel in Figure [5] shows the 
|ALPHA|-distribution of events with estimated energies be- 
tween 80 and 100 GeV; the right panel shows the distribution 
of events with reconstructed energies above 200 GeV. In both 
cases an excess of 7-ray events is clearly visible. However, 
the significance of the 7-ray signal at lower energies is con- 
siderably reduced compared to higher energies because of the 
degradation of the background suppression towards lower en- 
ergy (cf. Figure 0). 

At the lowest energies ALPHA is currently the only means 
by which it is possible to separate 7-rays and background 
events. This is illustrated in Figure|6] which shows the qual- 
ity factor as a function of SIZE separate for an optimized 
ALPHA cut and optimized HADRONNESS cut, as well as 
the combination of both cuts. Below SIZE 250 phe a cut in 



HADRONNESS does not improve 7/hadron-separation and 
reduces only statistics. On the other hand, with ALPHA a 
quality factor of 2 is still possible. 

4.3. Sensitivity 

The highest integral sensitivity for 7-ray emission from the 
Crab Nebula is obtained above ~ 250 GeV. The integral sensi- 
tivities for energies above 250 GeV obtained on a daily basis 
are listed in Table [2] The sensitivities are calculated in units 
of QLiMa /(hours) 1 / 2 ; the sign i ficanc e <7LiMa is calculated using 
equation 17 from lLi and Mai (119831) . In the last column of the 
table, the sensitivity is expressed as the minimum flux, nor- 
malized to the Crab flux that can be measured with 5 a signif- 
icance 24 in a 50 hr observation. The day-by-day sensitivities 
vary by about 10%, indicating a stable telescope performance 
throughout the observations. 

The energy dependence of the sensitivity was studied by 
calculating the integral sensitivity for several analysis thresh- 
olds. The SIZE-dependent HADRONNESS cut (Figure g] 
black solid line) was used for 7/hadron separation and the 
ALPHA cut was tightened with increasing energy based on 
studies from MC-simulations. The integral sensitivity for a 
50 hr observation is 13% Crab at 75 GeV and improves con- 
tinuously with increasing the analysis threshold to about 2.2% 
above ~ 250 GeV (see Table Note that above 1.6 TeV the 
background is estimated from only two events, which results 
in an uncertainty of more than 70% on the integral sensitivity. 

TABLE 2 

MAGIC Integral Sensitivities to 7-RAY Emission 
from the Crab Nebula for 7-RAY Energies above 250 GeV 



Date Sensitivity in Sensitivity in 

(MJD) [ouMa/Chours) 1 / 2 ] (% Crab, 50 hr, 5a) 

53,648 liU 2A 

53,655 19.6 2.1 

53.671 19.9 2.1 

53.672 19.1 2.2 

53,679 19.4 2.2 

53,684 21.1 1.8 

53,707 23.6 2.0 

53,709 16.1 2.7 

53,711 18.9 2.2 

53,713 18.7 2.2 

53,727 19.1 2.3 

53,729 17.5 2.3 

53,731 18.0 2.3 

53,735 18.7 2.1 

Average 19.1 ±0.5 (1.7) 2.21 ±0.05 (0.2) 

NOTE. — Average values of each column are given in the last row with 



the corresponding RMS value in brackets. Note that different definitions of 
significance are used in the two columns. 

4.4. Systematic Uncertainties 

Apart from statistical errors, many results in cosmic -ray ex- 
periments are affected by rather large systematic errors. One 
of the main problems is the lack of test beams that allow 
calibration of the entire instrument in combination with the 
showering process in the atmosphere. The standard replace- 
ment for a test beam calibration is the use of MC simulation 
based on many linked processes either using physical models 

24 Note that in this case the significance is calculated as 
excess / (background) '/ 2 . 
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FIG. 3. — Distribution of the image parameter WIDTH for 7-ray events extracted from data and for MC simulated 7-rays in four consecutive bins of SIZE 
covering the full analyzed energy range. 



(example: the simulation of electromagnetic showers) or cali- 
brating them in separate measurements (example: the spectral 
mirror reflectivity). In some cases one has to make reasonable 
guesses (example: photoelectron collection efficiency in the 
PMT front-end volume). The calibration of individual effects 
suffers partly from cross-correlations, which are not always 
well understood. Currently, the best approach is to estimate 
the systematic uncertainties (commonly called systematic er- 
rors) of the various parameters separately and to combine 
them to a global systematic error. Here, we follow the gen- 
eral practice of adding these individual errors in quadrature 
although this will result in a slight underestimate of the total 
systematic error. Table [9] lists the dominant systematic error 
contributions (> 2%) for the spectral parameters (flux, slopes, 
cut-offs, etc.). The systematic errors influence the spectrum in 
different ways. Some (Nr 1-7, 10, 11) result in an uncertainty 
of the energy scale and thus can enter as a large factor in the 



flux at a given energy in case of steep spectra with slopes < 3; 
others (Nr 8, 9, 12-13, 15-17) linearly influence the flux nor- 
malization for the spectrum. 

The most critical contributions to the systematic error come 
from the uncertainties in the conversion of photons to measur- 
able photoelectrons (combined under item 7), the so-called 
photon detection efficiency (PDE). The PDE is a combination 
of many small effects such as the reflectivity variation of the 
light catchers, tolerances in the light catcher geometry, angu- 
lar effects on the PMT surface, non-uniformity of the diffuse 
lacquer coating, the QE-spread and cathode non-uniformity 
of the PMTs, the photoelectron collection efficiency in the 
PMT front-end volume, and gain variations of the first dyn- 
ode. Also some contribution of the signal transmission to the 
DAQ is included. 

Fortunately, the PDE can normally be measured with a light 
source uniformly illuminating the camera with short blue or 
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FIG. 5. — Distribution of the image parameter |ALPHA| for the bin of reconstructed energy 80- lOOGeV in the left panel and for events with energies 
> 200 GeV in the right panel, (black: ON-data; red: OFF-data) 



UV light pulses. Obviously, the light pulser itself introduces 
some systematic errors such as in the absolute light flux de- 
termination, small deviations from uniformity in the illumi- 
nation, some (small) temperature drift, and amplitude jitter. 
Also, the used method of determining the number of de- 
tected photoelectrons, the above-mentioned so-called F-factor 
method, introduces some uncertainty. 

Another rather big uncertainty is the effective reflectivity 
of the mirrors defined as the light from a source at infinity 
being focused onto the area of a pixel. Comparing the mea- 



sured brightness of a star and its image back-reflected by a 
high-quality diffuse reflector in the camera plane allows one 
to carry out a routine reflectivity measurement (see Figure|7]), 
with an uncertainty of about 7%. A similarly large error con- 
tribution was estimated for the event reconstruction. Again, 
many small effects contribute to the reconstruction losses or 
to the wrong assignment of events. In contrast to the proce- 
dure to limit the uncertainty as in the example of the PDE, 
no simple method to cross-check the error range of the recon- 
struction efficiency is possible and a reasonable guess had to 
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TABLE 3 

Integral Sensitivities of the MAGIC Telescope to the 7-Ray Emission 
from the Crab Nebula for Several Analysis Thresholds 



Energy |ALPHA| Cut Sensitivity 

(GeV) (deg) ON Events OFF Events Excess Events (% Crab, 50 hrs, 5a) 

>75 < 10° 232505 221751 10754 131) 

> 110 < 10° 58600 49702 8898 7.5 

>200 < 10° 11399 4960 6439 3.3 

>400 <7.5° 2866 348 2518 2.2 

> 800 ■ 5 613 20 593 2.3 

> 1600 <5° 43 2 41 1.0 
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FIG. 6. — Size dependence of the Q-factor for optimized | ALPHA cut and 
HADRONNESS cut, and the Q-factor for the combination of ALPHA | cut 
and HADRONNESS cut. 

be made. 

Effects that influence the slope of reconstructed 7-ray en- 
ergy spectra (class C effects in Tabled are mostly dominating 
at the lowest and highest accessible energies. The estimate 
of the systematic slope error is rather difficult. In case of a 
power law or moderately curved power law we estimate an 
uncertainty on the slope of 0.2. We note that measurements 
by current second-generation telescopes of the spectral slope 
of the Crab Nebula agree better than 0.1 in the overlapping 
energy range. 

In summary, we obtain a systematic energy scale error of 
16%, a systematic error of 11% on the flux normalization 
(without the energy scale error), and a systematic slope er- 
ror of ±0.2 (which is a combination of error 13 and the other 
relevant class A errors averaged over the energy). 

5. ANALYSIS RESULTS 

5.1. Differential Energy Spectrum of the Crab Nebula 

By extracting the 7-ray signal in each bin of the recon- 
structed energy E rec , a spectrum Ni of 7-rays in each E lec bin 
i can be constructed. The reconstructed energy is subject to 
a bias. Before determining a differential 7-ray flux in true 
energy E lme bins, the spectrum N/ has, therefore, to be con- 
verted into a spectrum Mj of 7-rays in bins of Etrue- This is 
done by applying an u nfolding procedure with regularization 
dAnvkevev et alj|1991l) . An essential input for the unfolding 
procedure is the migration matrix, which describes the migra- 
tion of events from bin i in E KC into a bin j of £, rue . The mi- 
gration matrix is determined from MC simulated 7-ray show- 
ers. The unfoldin g is done independently for different regular- 
ization schemes (iTik honov and Arsenin ll979t lBerteroiri989l : 
ISchmelIing||1994t lAlbert et al.ll2007d) . Figure [8] shows one 
distribution of excess ev ents from the C rab Nebula after un- 
folding by the method of Bertero (1989). The integral rate of 



excess events is 0.4 Hz. The differences between the unfolded 
points Mj obtained with the different regularization schemes 
are used to estimate a systematic error due to the unfolding. 
Figure [10] shows the differential 7-ray flux, which was ob- 
tained with the regularization scheme proposed by iBerterol 
( 1989) and by normalizing the unfolded spectrum Mj, to the 
effective collection area (Figure [5), the effective observation 
time, and the bin width of Etnie (given by the horizontal bars 
at each flux point in the figure). The average differential flux 
for each energy bin is presented in Table [4] 

The influence of different choices for tail-cuts, HADRON- 
NESS cuts, DIST cuts, and core-pixel cuts on the measured 
flux is indicated in Figure [l0]by the shaded region and quoted 
as systematic uncertainty in Table |4] Due to analysis uncer- 
tainties, the band broadens at low energies, mostly because 
of limited 7/hadron discrimination power. It broadens at the 
highest energies due to low event statistics. 

The energy spectrum is parameterized with both a power- 
law and a curved power-law Ansatz. The fit takes into account 
correlations between the spectral points that are introduced by 
the unfolding procedure. A correlated fit with a power-law 

dF r 

— =/o(£/300GeV) , (1) 

provides a flux normalization ff of (5.7 ± 0.2 stat ) x 
lO-^cirrVTeV -1 and a spectral index T of -2.48 ± 
0.03 sta t ± 0.2 syst . The \ 2 of the fit is 24 for 8 degrees of 
freedom, which disfavors a pure power-law description of 
the spectrum. The energy spectrum is better described by a 
curved power-law Ansatz 

g = ff (£/300GeV) [fl+M ° g >° (£/30 ° GeV) ] (2) 

yielding a flux normalization ff of (6.0 ± 0.2 stat ) x 
lO-^cnrV'TeV- 1 , a = -2.31 ± 0.06 stat and b = -0.26 ± 
0.07 stat ±0.2 syst . The \ 2 of the fit is 8 for 7 degrees of free- 
dom. 

Figure [TTJ shows the differential flux measurements multi- 
plied by the energy squared, i.e. the spectral energy distribu- 
tion (SED). In the figure we compare our measurement with 
those from other experiments. For energies above 400 GeV 
the derived spectrum is in good agreem ent with measure- 
ments of other air Cerenkov telescopes ([Hillas et al.l Il998h 
l Aharonian et al.ll2004t iTanimori et alJl 1998b lAharonian et al] 
120061) . At energies < 400 GeV, below the threshold of 
previous measurements by IACTs, we compare our results 
with those obtain ed by CELEST E dde Naurois etail 120021: 
ISmith et alJl2006l) and STACEE (lOser et al.ll2001l) . i.e. mea- 
surements performed by converted solar tower experiments. 
It should be noted that the integral flux values of these exper- 
iments had to be converted to differential ones by assuming 
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FIG. 7. — Measurement of the reflectiv ity of the main mirr or; the effective 
mirror area is shown by the blue points I Garczarczyk 2006). The measure- 
ments cover the Crab observation period from 2005 October to December. 
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FIG. 9. — Collection area after image cleaning and after cuts for low zenith 
angle observations (< 20°). 



a shape of the source spectrum, which causes an additional 
bias. 

Above 1 TeV up to about 10 TeV the measure d energy spec- 
trum is well described by a pure power-law dWeekes et al.l 
Il989t lAharonian et al.ll2004l) . Going to lower energies, one 
expects a continuous hardening of the spectrum. However, 
this could not be demonstrated by earlier measurements. The 
change of the slope of the spectrum (spectral index) r" was 
tested for various points of the measured spectrum r". 

dln(F) Aln(F) InFj-lnFj 

T'(E)= — w- — — w ■ J - (3) 

dln(E) Aln(E) ln£,-ln£ ; v ; 

E = exp [0.5(ln£; + ln£;)] , (4) 

where Fy is the differential flux measured at Eij. The four 
derived spectral indices at ~ 150 GeV, ~ 300 GeV, ~ ITeV, 
and ~ 2.5 TeV shown in Figure Q~2] indie ate a clear softening 
of the spectrum with increasing energy. The spectral index 
T' was also derived from the aforementioned results of the 
curved power-law fit, 

r' = a + 2Mog 10 (£/300GeV) (5) 

and is shown by the solid black line, and the ± 1 er-confidence 
band is shown by the dashed black line. A systematic uncer- 
tainty on the slope can cause an additional vertical shift of the 
measurement by ±0.2. Within uncertainties, the measured 
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FIG. 8. — Unfolded distribution of excess events from the Crab Nebula. 
The integral rate of excess events is 0.4 Hz. 



spectral index varies in good agreement with predictions by 
lAharonian et al.l d2004l) (blue line), who, in addition to the IC- 
scattering on synchrotron photons, included several other soft 
photon fields such as millimeter-photons, CMB, and far-IR 
photons from dust and stars. 

The predicted GeV 7-ray emission has a peak in the SED- 
representation (see Figure fTTTi. If one assumes that the en- 
ergy spectrum around the peak can be described with a curved 
power-law, the position of the peak can be determined from 
the measurement of the spectral index obtained from the re- 
sult of the curved power-law fit. A necessary condition for 
the peak in the SED is that the spectral index T' is -2. With 
this condition the peak is determined at 77 ± 47 stat "^J ^ GeV 
(Figure [12] triangle). 

5.2. VHE ry-Ray Light Curve of the Nebula Emission 

In the VHE 7-ray astronomy community it is assumed that 
the Crab Nebula is a constant and stable 7-ray source and can 
therefore be used as a standard candle. However, with more 
sensitive measurements it is necessary to check the stability of 
the 7-ray source. Below we present a time-resolved measure- 
ment of the VHE 7-ray flux, i.e. the light curve for the Crab 
Nebula. Depending on the source strength and the analysis 

TABLE 4 

Mean Energy and Differential Flux of the Spectral Points 
Shown in Figure[To] 



EnergyMean Energy 


Differential Flux 


(GeV) 


(cnrV'TeV- 1 ) 


77 


(1.14±0.27 sM ±0.34 syst )10- 8 


127 


(3.65±0.38 s tat ±0.55 sys t)10- 9 


210 


(1.41±0.09 slal ±0.28 sys ,)10- 9 


346 


(4.37±0.23 slal ±0.87 syst )10- 10 


570 


(1.32 ± 0.07 stal ± 0.20 syst ) 10-'° 


940 


(3.55 ± 0.23 slal ± 0. 18 syst )10-' 1 


1550 


(9.88±0.74 slal ±0.49 syst )10- 12 


2554 


(2.69±0.29 slal ±0.27 syst )10- 12 


4212 


(6.80 ± l.lOstat ± 1.00 sys ,)10- 13 


6943 


(1.15±0.53 sla , ±0.12 sys ,)10- 13 



NOTE. — The systematic errors are derived from different applied cuts 
and unfolding procedures. 
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was unfolded with the method of Bertero 1 1989). The results of a fit of the 
spectrum with a power law and a broken power law are also shown. The 
bottom panel shows the relative residuals between the fit and the data points. 
See text for further discussion. 

threshold, the time intervals can be as short as a few minutes. 

We calculated light curves in bins of 10 minutes from events 
with estimated energies above 200 GeV. The light curves of all 
14 selected nights are shown in FigureQj] Note that the same 
loose cuts are used for the 7/hadron-separation as for the cal- 
culation of the energy spectrum, which reduces the sensitivity 
of the measurement. The probability that the light curve is 
described by a constant flux level is > 10% in all nights ex- 
cept the first night, where the probability of the fit is 0.8%. 
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FIG. 11. — SED of the 7-ray emission of Crab Nebula. The measure- 
ments shown below 10 GeV are by EGRET jde Jager et aTiri996l) . In VHE 
7-rays, measurements are from ground-based experiments. Above 400 GeV 
our measurement is in agreement wi th measurement s by other IACTs. The 
dashed line is a model prediction by Aharonian et al. 12004). 



MAGIC, derived from flux points 
MAGIC, derived from variable power law fit 
68% confidence interval 
MAGIC, peak position 
Aharonian et al. (2004) 
Atoyan & Aharonian (1996) 
Atoyan & Aharonian (1998) 




-3.4 



10 2 1 1 4 1 

Energy [GeV] 

FIG. 12. — Measured spectral index derived from differential flux 
points (filled circles) and from the curved power-law fit (black solid 
line; the dashed line gives the ltr confidence interval); Predictions by 
lAharonian et al.l 420041) (blue line), Atoyan and Aharonian (1996) (green 
line), and Aharonian and Atovan 1 1998) (red line) are also shown. 

The average statistical uncertainty of each flux measurement 
is ~ 20%. Figure [14] shows the average flux of each night. 
The dashed line in the figure denotes the average flux from all 
nights, and the shaded region shows the statistical error in the 
flux. The average integral mean flux f>200GeV is 



f>200Gev = (1.96±0.05 smt )xlO- 



10 cm- 2 s- 



(6) 



There is a probability of 67% that the measured daily flux 
values are compatible with a constant flux. We can, there- 
fore, conclude that the reconstructed flux of the Crab Nebula, 
within statistical uncertainties, was constant over the entire 
observation period. 

5.3. Morphology of the ^y-Ray emitting Region 

The morphology of the 7-ray emission was studied by gen- 
erating sky-maps in three uncorrelated bins of SIZE. The 
reconstruction of the origin of a 7-ray event with a sin- 
gle telescope is possible with the so -called DISP-method 
(iFomin et al.lll994t iLessard et alJl2001|) . For the studies pre- 
sented here we used the following parameterization for DISP: 

WIDTH 

DISP = a (SIZE) + b (SIZE) , (7) 

V ' V ; LENGTH ' V ; 

where a and b are second-order polynomials found by fit- 
ting M C simulated 7-ray showers dDomingo-Santamarfa et alj 
2005). Strong tail cuts of 10 and 8 phe were used in the im- 
age cleaning for core- and boundary-pixels, respectively, and 
a tight HADRONNESS cut < 0.1 was applied, resulting in 
improved angular resolution. 

The reconstructed event origins were corrected for possi- 
ble mispointing by using the information from the starguider 
camera. Two-dimensional (2D) histograms with bin sizes of 
(0.057° x 0.057°) were filled with the corrected event origins 
(events with energies < 500 GeV). A 4 times finer binning 
was chosen for the sky-map filled by events with energies 
above 500 GeV. Figure [151 shows the background-subtracted 
sky-maps of excess events from the Crab Nebula for 7-ray 
energies - 160, - 250, and > 500 GeV. 

5.3.1. Center of Gravity of the -y-Ray Emission 
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FIG. 13. — Light curves of the integral Flux above 200 GeV from the Crab Nebula for each night. 



The CoG of the 7-ray emission was derived from the sky- 
maps of the excess events shown in Figure[l5]by fitting them 
with a 2D-Gaussian of the form 

(x-x) 2 + (y-y) 2 



F(x,y) = F les + a-e\p 



2a 2 



(8) 



where f les is introduced to account for a possible constant off- 
set of the background-subtracted sky-map. In this representa- 
tion a defines the 39% containment radius of the observed 7- 
ray emission. Here we assume that the distribution of excess 



events is rotationally symmetric, i.e. a x = a v = a. It is further 
assumed that a is the convolution of the response of the de- 
tector <7 ps f (point-spread function) and the apparent size of the 
7-ray emission region a sm , i.e. a 2 = a 2 f + of rc . Note that here 
it is assumed that the 7-ray emission region follows the shape 
described by Equation [8] which in reality is not necessarily 
the case. 

The CoGs obtained from the fitted x and y are listed in Ta- 
bleland shown in Figure [16] superimposed on the composite 
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MJD 
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53731.0 
53729.0 
53727.0 
53713.1 
53711.1 
53709.1 
53707.1 
53684.1 
53679.2 |- 
53672.1 
53671.2 
53655.2 
53648.2 
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FIG. 14. — Nightly average flux from the Crab Nebula above 200GeV of 
each observed night. The dashed blue line gives the average flux of all nights 
and the blue shaded region gives the corresponding statistical error. 

image of optical, IR, and X-ray observations of the Crab Neb- 
ula. The three measured CoGs are compatible among each 
other and coincide with the position of the pulsar. Note that 
the systematic uncertainty of the position is ~ 1'. 

5.3.2. Extension of the j-Ray Emission Region 

The extension of the 7-ray emission region was studied by 
comparing the width of the excess event distribution with that 
obtained for a simulated 7-ray point-source. The simulated 
distributions were verified by comparing them to the distribu- 
tions extracted from an observation of Mrk 42 1 , an extragalac- 
tic 7-ray source that can be considered a point so urce for our 
purpose . The Mrk 421 data set is the same as in Alb ert et all 
(2007a). The width of the 7-ray excess extracted from Mrk 
421 and the simulated width for a point source agree within 
statistical uncertainties. 

In the following, the average position of the CoGs obtained 
from the three sky-maps is assumed as the 7-ray source po- 
sition. The angular distance squared (8 2 ) between the recon- 



structed origin and the assumed source position is calculated 
for every event. The background subtracted (^-distributions 
obtained for the three energy ranges are shown in Figure [TT] 
Data (black) and MC (blue) are compatible within statistical 
uncertainties in all three (^-distributions. 
An exponential function of the form 



(9) 



describes the expected f? 2 -distribution, where a is a normal- 
ization and g is the same as in equation[8] Values for a 1 and 
fip sf found by fitting the corresponding 8 2 -distributions with 
equation [9] are shown in Table [6] Upper limits on (j src were 
calculated with a confidence level o f 95% following the pro- 
cedure outlined in I Yao et"ail (120061) for one-sided confidence 
intervals and Gaussian errors. The results are presented in Ta- 
ble[6j For energies above 500 Ge V the limit is shown in Figure 
[161 The limits obtained for 7-ray energies above 500 GeV and 
about 250 GeV constrain the 7-ray emission to a region within 
the optical synchrotron nebula. 

5.4. Search for pulsed ^-Ray Emission 

Among the most challenging tasks of ground-based 7-ray 
experiments is the detection of a pulsar. Several experiments 
have tried but failed. Currently MAGIC is the only ground- 
based detector with threshold settings below 100 GeV that is 
appropriate for a search of pulsed 7-ray emission from the 
Crab pulsar. For the data a periodicity analysis was performed 
after 7/hadron-separation and selection of events with small 
I ALPHA I -value. The cuts were chosen by MC simulations to 
optimize the sensitivity of the analysis. After event selection, 
the event times 25 were transformed to the baryc enter of the 
solar system with the TEMPO timing package (Tayl or et alj 
2000, 26 ). Then, the corresponding phase <pj of the Crab pul- 
sar was calculated for each transformed arrival time t ,•: 



1 ■> 
b j = v(tj-k ) ) + -v(t j -t Q ) 2 



(10) 



where v, v, and fo are values of contemporary ephemerides 
of the Crab pulsar provided by the Jodrell Bank Radio Tele- 
- 17 (see Table IT). W e tested for periodicity with the H- 



scope 

test (Ide Jager et al 



1989), the Pearson's x 2 -test, and a test 



TABLE 5 

Center of Gravity of the 7-RAY Emission of the Crab Nebula 
Obtained for Different Energies 



SIZE 
(phe) 


Energy 
(GeV) 


Right Ascension 
(hr) 


Declination 
(deg) 


200-300 
300-700 
> 700 
Position of the pulsar: 


160+ 80 
> 500 


5.5766 ±0.0009 

5.5758 ± 0.0003 

5.5759 ± 0.0003 
5.5755 


22.019 ±0.011 
22.019 ±0.004 
22.022 ±0.003 
22.015 



NOTE. — In the first column the applied SIZE cut is stated. The second 
column shows the corresponding range of 7-ray energies covered (peak value 
and full width at half maximum of the distribution of MC 7-ray events). The 
last two columns give the fitted position of the CoG and the statistical uncer- 
tainty. 

25 The time of each event was derived from the time signal of a GPS con- 
trolled rubidium clock with a precision of ~ 200 ns. 

26 Available at http://www.atnf.csiro.au/research/pulsar/tempo/ 



See http://www.jb.mac.ac.uk/pulsar/crab.html 
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FIG. 15. — Sky-maps of excess events from the Crab Nebula for different 7-ray energies (~ 160, ~ 250, and > 500 GeV). The position of the pulsar is marked 
by the black plus sign and the angular resolution is indicated by the circle. 



from iGregorv and Loredol (119921) that is based on Bayesian 
statistics. 

The analysis chain was tested by optical observations of the 
Crab pulsar with the MAGIC telescope. Within this 12.5hr 
observation, every time the readout of MAGIC was triggered 
by a cosmic-ray shower, the signal of the pixel in the center 
of the MAGIC camera was recorded by the MAGIC DAQ for 
100 ns. Along with an average trigger rate of 200 Hz, the ef- 
fective observation time was only about 1 s. Figure [T8l shows 
the reconstructed optical light curve of the Crab pulsar with 
the familiar main pulse and interpulse. For better readability 
the light curve is shown twice. The position of the main-pulse 
is shifted with respect to the position of the main-pulse in ra- 
dio by -252 ± 64/is, w hich is in agreement with the contem- 
porary measurement of Oosterbroek et al. (2006). 

5.4. 1 . Search for pulsed Emission in Differential Bins of Energy 



We searched for pulsed 7-ray emission in five bins of re- 
constructed energy between 60 GeV and 9TeV. This search 
was motiv ated by a possible p ulsed 7-ray component at TeV 
energies (Hirotani 2001, 2007). However, no signature of pe- 
riodicity was found in any of the tested energy intervals. 

For each energy bin an upper limit on the number of ex- 



TABLE 6 

Results of the fit of the ^-distributions with an exponential 
ansatz and thereof derived upper limits on the extension of 
the emission region (39% containment radius). 



Energy 


Data <j~ 




95% U.L. on Extension 


(GeV) 


(deg 2 ) 


(deg 2 ) 


(arcmin) 


160+so 


0.0148 ±0.0035 


0.01 13 ±0.0007 


5.9 


250 +i30 


0.0100 ±0.0008 


0.0100 ±0.0003 


2.4 


> 500 


0.0054 ±0.0006 


0.0051 ±0.0002 


2.2 
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TABLE 7 

Ephemerides of the Crab Pulsar from the Jodrell Bank 
Telescope Covering the Same Period of Time as the Analyzed 

Data 



t Frequency (u) Derivative (u) 

JD (s) (Hz) (1(T 15 Hzs" 1 ) 

2,453,597.5 0.029626 29.7798524524 -372992.36 

2,453,628.5 0.031767 29.7788534525 -372972.07 

2,453,658.5 0.022656 29.7778867428 -372950.45 

2,453,689.5 0.016803 29.7768878849 -372924.54 

2,453,719.5 0.026788 29.7759213143 -372886.52 

2,453,750.5 0.020341 29.7749226318 -372854.62 

2,453.781.5 0.006520 29.7739240139 -372823.70 

NOTE. — Given in JD is the reference day of the ephemeris, and ; is the 
time of appearance of the first main pulse on the reference day after midnight. 




FIG. 16. — Emission of the Crab Nebula in different energy bands. The 
position of the Crab pulsar is marked with a black star. The Chandra X-ray 
image is shown in light blue, the Hubble Space Telescope optical images are 
in green and dark blue, and the Spitzer Space Telescope's infrared image is in 
red (picture from ?). Overlaid are the CoG of the 7-ray emission at different 
energies (plus sign: > 500 GeV; filled circle: ~ 250 GeV; filled square: ~ 
160 GeV). The error bars indicate the statistical uncertainty in the position of 
the CoG. Indicated by the dashed circle is the upper limit (95% confidence 
level) on the 39% containment radius of the 7-ray emission region that was 
derived from the 9 2 -distribution for 7-ray energies above 500 GeV. 



cess events was calculated with a confidence level of 95% in 
two differ ent ways, first fr om the result of the //-test as de- 
scribed by lde Jag er ( 1994), and second from the pulse phase 
profile. In the calculation of the limit from the result of the 
//-test, it is assumed that the duty cycle (FWHM) of the 
pulsed 7-emission is 20%, similar to the duty cycle of the 
light curve measur ed by the EGRET detector above lOOMeV 
dFierro et al.lll998f) . No assumption about the position of the 
emission in the pulse phase profile enters into the calculation. 

This additional constraint is applied, however, when the up- 
per limit is directly derived from the pulse phase profile. As 
signal regions we chose the phase intervals where EGRET had 
observed pulsed emission above 100 Me V, i.e. -0.06 — 0.04 
and 0.32 — 0.43 (shaded region in Figure |20] below). The 
background was estimated from the remaining phase inter- 
vals. Having defined the signal and background regions in 
this way, the upper limit on the number of ex cess events was 
obtained by the method of lRolke et alj (|2005). Because of the 
additional constraint made about the position of the expected 
pulsed emission, the limits obtained from the pulse phase pro- 
file are on average about a factor of 2 better than the limits 
obtained from the result of the //-test. 

The upper limits derived from the pulse phase profiles were 
converted into flux limits. The collection area was calculated 
assuming a photon index of -2.6 for the 7-ray spectrum. The 
flux limits are shown in Figure [19] together with the upper 
limit on the cutoff energy, which is derived in the following 
section. 

5.4.2. Upper Limit on the Cutoff Energy of the pulsed Emission 

Apart from the search for pulsed emission in bins of re- 
constructed energy, we performed a periodicity analysis, this 
time selecting events with SIZE< 300phe (7-ray energies < 
180 GeV) and applying the same optimized SIZE-dependent 
HADRONNESS cuts and ALPHA cuts as above. Compared 



to the previously described analysis, this one is optimized for 
a search of pulsed emission close to the threshold of the ex- 
periment. The analysis threshold, defined as the peak of the 
energy-distribution of simulated 7-ray showers, is 60GeV. 

Figure [20] shows the pulse phase profile obtained for 
the selected events. For comparison the pulse phase pro- 
files from EGRE T observations above 100 MeV and 5 GeV 
(iThompson et al.l [20041) are al so shown. The EGRE T pulse 
phase profile above 10 GeV ( IThompson et al.l 120051) is not 
shown because it suffers from too low statistics. Shaded in 
the pulse phase profiles are the regions of the main and in- 
terpulse defined fr om EGRET observations above 100 MeV 
dNolanet alJll993h . 

The result of a Pearson's x 2 -test is 13.1 with 10 degrees of 
freedom, corresponding to a significance of 1 .2 a for periodic 
emission. The result of the //-test is 3.9 , which is equivalent 
to a si gnificance of 1.3 a. The test by iGregorv and Loredol 
(1992) results in a probability of 4.1 ■ 10 4 that pulsed emis- 
sion is present in the data. These tests do not make an assump- 
tion about the position of the pulsed emission in the pulse 
phase profile. However, some evidence of an excess is visi- 
ble at the position of the inter-pulse in the same phase range 
where EGRET detected pulsed emission above 100 MeV. If 
the two phase regions defined by EGRET are used as the sig- 
nal region and the remaining phase intervals as background 
region, the significance of the observed excess is 2.9 a. Note 
that in this case the significance was not calculated from the 
binned pulse phase profile shown in Figure |20l 

The significance of the observed excess is not sufficient to 
claim the detection of a pulsed signal; therefore, upper limits 
on the number of excess events were calculated with a con- 
fidence level of 95% (see Tabled. Note that because of the 
observed excess, the upper limit from the pulse phase pro- 
file is larger than the limit obtained from the //-test. Using 
the different limits on the number of pulsed excess events, 
we constrain, in the following, the cutoff energy of the pulsar 
spectrum under the assumption that the break in the energy 
spectrum can be described with an exponential cutoff. In the 
procedure we use the par ametrization of th e measured pulsar 
spectrum below 10 GeV dFierro et alJI 19981) . extended with an 
exponential cutoff: 



F(E , -Ecutoff ) : 
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FIG. 17. — Background subtracted 9 2 -distributions for different energies. ~ 160GeV (top left), ~ 250GeV (top right), > 500GeV (bottom) 



x exp 



^Cutoff / 



photons (cm 2 sGeV) (11) 



in Tabled] 



The spectrum with a given £cutoff is convoluted with the 
effective collection area after cuts. The collection area is de- 
rived from MC simulations, assuming the same 7-ray spec- 
trum. The number of expected excess events for the assumed 
cutoff energy is obtained by multiplying the convoluted spec- 
trum with the observation time. In an iterative algorithm 
^Cutoff is changed until the number of expected excess events 
matches the upper limit on the number of excess events. In 
this way we derive an upper limit on the cutoff energy of 
27GeV from the result of the //-test and 34GeV from the 
limit obtained from the pulse phase profile. Differential and 
integral upper limits on the flux were calculated and are shown 
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FIG. 18. — Optical light curve of the Crab pulsar measured with MAGIC. 
The figure includes data from seven different observations between 2005 De- 
cember and 2006 February. The total observation time was 12.5 hr. 



6. DISCUSSION 



In this paper we report on the most detailed study to date of 
VHE 7-ray emission of the Crab Nebula below 500 GeV. This 
study includes the following: 

• A measurement of the differential energy spectrum 
down to 60 GeV. 
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FIG. 19. — Upper limits on the pulsed 7-ray flux from the Crab pulsar; 
upper limits in differential bins of energy are given by the blue points. The 
upper limit on the cutoff energy of the pulsed emission is indicated by the 
dashed line. The analysis threshold to derive the upper limit on the cutoff 
energy is indicated by the red arrow. 
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TABLE 8 

Analysis Results for a Cut selecting Events with SIZE < 300phe 



Method //-Test Rolke 



Test result 3.92 n.a. 

Significance 1.26<r n.a. 

2<r U. L. on excess events 1635 3198 

U. L. on the Cutoff energy (GeV) 27 34 

2<r Integral flux limit above 60 GeV (cnr 2 s _1 ) 2.5 ■ KT 1 1 7.9 ■ 10"' 1 

2a Differential flux limit at 60 GeV (cirT 2 sT 1 GeV" 1 ) 4.5 • 10~ 12 8.9 ■ 10~ 12 

Peak Energy MC (GeV) 60 60 



• An estimate of the peak in the SED of the VHE 7-ray 
emission. 

• The search for an extended source morphology. 

• The calculation of a light curve of the VHE 7-ray emis- 
sion from the nebula above 200 GeV. 
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FIG. 20. — Pulse phase profiles of the Crab pulsar. Bottom panel: opti- 
cal observations by MAGIC (~ 1 s effective observation t ime); middle pan- 
els: observations by EGRET from Thompson et al. (2004); top panel: pulse 
phase profile obtained by MAGIC, analysis threshold 60 GeV. The shaded re- 
gions indicate the EGRET measured positions of the peaks for 7-ray energies 
above lOOMeV. Note that each pulse phase profile is shown twice for better 
visibility. 



• A search for pulsed emission from the Crab pulsar in 
differential bins of energy and in an optimized low- 
energy analysis. 

Most of the aforementioned studies were done in this energy 
region for the first time; they were possible only because the 
imaging air shower Cerenkov technique was used. The wave 
front sampling technique, which, up to now, was the only 
experimental technique used in this energy domain, did, at 
best, allow one to arrange an integral flux measurement and 
to search for pulsed emission. The performance of MAGIC 
is superior even in the energy domain below 200 GeV, where 
a progressive degradation of the 7/hadron-separation power 
is observed. In this context studies for improving the sup- 
pression of background events by exploiting the intrinsic time 
structure of the recorded PMT-signals are ongoing. Further 
improvement is expected with the second MAGIC telescope 
currently under construction. 

The measured energy spectrum of the Crab Nebula in Fig- 
ure[lO]extends over two decades in energy and five decades in 
flux. The spectral shape deviates from a pure power-law be- 
havior and is, within experimental uncertainties, in agreement 
with a curved power-law. The observation supports the gener- 
ally accepted picture that the steady emission above a few tens 
of GeV and up to the highest measured energi es can be de- 
scribed within the framework of the SSC-model (iGouldl 1965b 
Ide Jager and Harding|[T99l lAharonian et al.ll2004l) . Also the 
peak position of the inverse Compton emission of the tested 
predictions is in agreement with the estimated peak in the 
SED (77 ±47 stat + _^ 7 syst GeV). 

At GeV energies EGRET observed a 7-ray flux, which was 
a factor of 5 above th e flux predicted by the S SC-mechanism 
dde Jager et al.lll996l) . lAtovan and Aharonianl d 19961) explain 
this 7-ray excess by an additional 7-ray component from 
bremsstrahlung of electrons that are partially captured in fil- 
aments of the nebula. Such an extra component can signifi- 
cantly change the spectral slope at several hundred GeV com- 
pared to a pure IC-scenario (cf. blue and red lines in Figure 
[T2l and results in an almost pure power-law behavior of the 
energy spectrum between ~ 100 GeV and lOTeV (constant 
T). At several hundred GeV, where the measurement is most 
sensitive, the measured slope {black line and data point) is 
consid erably harder than predicted by Atova rTand Aharonianl 
( 1996). It is, therefore, unlikely that the 7-ray excess at GeV 
energies can be exp lained by bremsstrahlung as proposed. 
Later predictions by lAharonian and Atoyanl (119981) that also 
include the mentioned bremsstrahlung mechanism, are in 
agreement with the presented measurement (red line). How- 
ever, all above-mentioned predictions agree with the mea- 
surement if the measured slope is shifted by 0.2 to more 
negative values, which is within the range of the system- 
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atic uncertainty of the measu rement. In the prediction by 
lAharonian and Atoyanl 0998) also a 7-ray component from 
7r°-decay is included, which results in a considerable harder 
spectrum above a few TeV compared to the pure IC-scenario 
(cf. Figure [121). However, given the limited statistics above 
1 TeV of our measurement, one cannot exclude any such pre- 
diction from the measurement. 

Studies about the morphology of the 7-ra y-emitting region 
of the Crab Nebula have been performed by lAharonian et al.l 
(2000, HH for 7-ray energies above 1 TeV. In both cases 
it was found that within the resolution of the experiment the 
emission region is pointlike. They placed an upper limit on 
the source size of ~ 2' at energies between 1 and 10 TeV. In 
the VHE domain, the morphology of the emission region has 
not yet been studied at energies below 1 TeV. With the resolu- 
tion of MAGIC it was possible to constrain the origin of the 7- 
ray emission to be within the optical synchrotron nebula (see 
FigurefTo*]). The upper limit on the size of the emission region 
is ~ 2', which is about 4 times larger than the predicted size 
of the inverse C ompton surface brightness f or 7-ray energies 
below 500 GeV dde Jager and Hardinglll992l) . 

X-ray observations indicate variabilities in the acceleration 
and cooling times of electrons on timescale of months (e.g. 
iHester et al. 2002). However, variations in 7-rays could not 
be detected so far. The sensitivity of MAGIC allowed us to 
study the variability of the 7-ray emission above 150 GeV on 
timescales as short as a few minutes up to months. We mea- 
sured a flux that is within statistics compatible with steady 
emission. During the observation the stability of the integral 
flux was better than 10% on all tested timescales. 

In a search for pulsed VHE 7-ray emission with MC op- 
timized cuts in HADRONNESS and ALPHA an excess was 
found in the pulse phase profile at the same position where 
EGRET detected pulsed emission above 100 Me V. The sig- 
nificance of the excess is 2.9 a if the phase regions where 
EGRET detected pulsed emission were chosen as signal re- 
gions and the remaining phase intervals are considered as 
background regions. The similarity of the distribution of ex- 
cess events in the EGRET > 5 GeV and MAGIC data and the 
monotonic increase of the number of excess events with in- 
creasing upper SIZE cut are strong indications that the ob- 



served excess is not a random fluctuation. 

With the result of the //-test an upper limit on the cutoff 
energy of 27 GeV was derived, assuming that the power-law 
spectrum of the pulsar at GeV energies is attenuated by an 
exponential cutoff. However, if the cutoff of the spectrum has 
a super-exponential shape, a cutoff energy almost as high as 
the analysis threshold (~ 60 GeV) cannot be excluded. 

With the derived upper limit we constrain not only the 7- 
ray emission from within the light cylinder but also the pre- 
dicted pulsed 7-ray emission in the unshocked wind region 
(Bogovalov and Aharonian 2000). The predicted 7-ray flux 
strongly depends (1) on the distance from the light cylinder 
where the kinetic energy-dominated wind forms and (2) on 
the wind's Lorentz factor. By comparing our observational 
limits on the pulsed emission w ith the predicted spectra by 
iBogovalov and A haronian (2000|), we can exclude the forma- 
tion of a particle-dominated wind within a few light cylinder 
radii. Following the argumentation of the same authors, the 
particle-dominated pulsar wind must therefore be formed far- 
ther out, most likely at distances of more than 30 light cylinder 
radii. 

Also, no pulsed emission was detected for energies above 
100 GeV, which could have its origin in IC-upscattering of IR 
photons within the l ight cylinder. Despite earl ier claims of 
a st rong componen t dHirotani and S hibata 2001), latest mod- 
els ((Hirotani 2007j) seem to disfavor a pulsed TeV compo- 
nent from the Crab pulsar due to dominant 7-7 absorption 
processes. In the future, detailed spectroscopic studies of 
the pulsed emission by, e.g., GLAST and ground-based ex- 
periments with lower thresholds and higher sensitivities like 
MAGIC II (under construction) and CTA (projected) will 
hopefully resolve the long-standing question of the origin of 
the pulsed emission. 
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TABLE 9 

Contribution to the Systematic Uncertainties 



Uncertainty 

Item Source of Uncertainty Class (%) Comments 



1 Parametrization of Atmosphere in MC-simulation A 3 Deviations due to yearly and daily pressure changes, deviations of real density distribution 

and standard atmosphere model 

2 Atmospheric transmission losses due to Mie scattering A, (C) 5 Lack of good measurements; short term unpredictable changes possible 

3 Incorrect NSB simulation A 3 MC assumes uniform NSB. Variations due to source location, air glow, variations 

due to manmade light. Stars in the FoV. 

4 Reflectivity of main mirror A 7 From measurements of reflected star images 

5 Variation of the useful mirror area A 3 Malfunctions of active mirror control resulting in focussing losses 

6 Day to day reflectivity changes A 2 Due to dust deposit variations and occasional dew deposit 

7 Photon detection efficiency of the PMT/lightcatcher system A, C 10-12 See text 

8 Unusable camera channels B 3 Dead PMTs (5-10 channels), problems in calibration (5-10 channels) 

9 Trigger inefficiencies B, C 4 Due to discriminator dead-time, baseline shifts/drifts, level differences trigger branch and FADC branch etc. 

10 Signal drift in camera due to temperature drifts A, C 2 Combination of PMT QE change (small), amplifier and optical transmitter drifts 

11 Camera flatfielding A,B 2 Calibration problem 

12 Signal extractor B 5 Complex effect due to trigger jitter (early pulses from PEs generated on 1st dynode) etc.; 

baseline jitter, shifts in FADCs 

13 cuts and methods used in the analysis B,C 5-30 Energy dependent, see discussion of differential energy spectrum 

14 Losses of events during reconstruction B(A) 8 Simplifications in MC simulation 

15 Estimate of BG under source B(A) 4 Camera nonuniformity not included in MC. Hadronic events not perfectly simulated in MC. 

16 Small tracking instabilities B 2 Source jitters around nominal camera position due to small tracking errors, small camera oscillations 

due to gusts etc, resulting in a wider signal spread than predicted by MC 

17 Nonlinearities in the analog signal chain (PMT-FADC C(A) 3-10 Saturation and nonlinearities of electronic and opto-electronic components 

NOTE. — Class A: contributions to the uncertainty on the energy scale. Class B: contributions to the uncertainty in the event rate. Class B(A): error contributes more to the leading term. Some of the uncertainties 



are energy dependent and are averaged. Class C: contribution affecting the spectral slope. 



